fs = 1000;        % 采样率 (Hz)
t = 0:1/fs:1;     % 时间向量，1秒
f = 5;            % 方波频率 (Hz)
x = square(2*pi*f*t);   % 生成方波（幅值 -1 和 +1）

plot(t, x);
xlabel('时间 (s)');
ylabel('幅值');
title('方波信号');
